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We discuss a three-flavor Nambu-Jona-Lasinio (NJL) type quantum field theoretical ap- 
proach to the quark matter equation of state (EoS) with scalar diquark condensate, isoscalar 
vector mean field and Kobayashi-Maskawa-'t Hooft (KMT) determinant interaction. While 
often the diquark and vector meson couplings are considered as free parameters, we will 
fix them here to their values according to the Fierz transformation of a one-gluon exchange 
interaction. In order to estimate the effect of a possible change in the vacuum pressure of the 
(N ■ gluon sector at finite baryon density we exploit a recent modification of the Polyakov-loop 

Q |i NJL (mPNJL) model which introduces a parametric density dependence of the Polyakov- 

' loop potential also at T = 0, thus being relevant for compact star physics. We use a 

TJ*± ' Dirac-Brueckner-Hartree-Fock (DBHF) EoS for hadronic matter phase and discuss results 

for mass-radius relationships following from a solution of the TOV equations for such a 
hybrid EoS in the context of observational constraints from selected objects. 

§1. Introduction 



In the present contribution we will introduce a model for the EoS of hybrid stars 
with deconfined quark matter cores and recent observational constraints^'^ for their 
masses (M) and radii (R). 

The status of the theoretical approach to the neutron star matter equation of 
state is very different from that for the high-temperature case at low or vanishing 
net baryon densities, where ab initio lattice QCD simulations provide EoS with al- 
most physical quark masses systematically approaching the continuum limit. 3 ' This 
guidance is yet absent at zero temperature and high baryon number densities, where 
a variety of models exists on different levels of the microphysical description which 
make different predictions for the state of matter. It is the current hope that the 
situation might change in a not too far future when very accurate measurements of 
the M — R relationship for compact stars become possible, e.g., with the Interna- 
tional X-Ray Observatory (IXO) project.® Once we are in possession of such data 
which allow upon inversion of the Tolman-Oppenheimer-Volkoff (TOV) equations 
the extraction of the cold dense compact star EoSP we will face the problem of 
interpreting the physical content of this numerical result. Then, a broad basis of 
alternative EoS models like the one presented here may become very useful. 

A common feature of present hybrid star models is that the transition from 
hadronic to quark matter cannot yet be described on a unique footing where hadrons 
would appear as bound states (clusters) of quarks and their possible dissociation at 
high densities as a kind of Mott transition® like in nonideal plasmas^ 1 or in nu- 
clear matter.^® Early nonrelativistic potential model approaches''^ 1 are insuffi- 
cient since they cannot accomodate the chiral symmetry restoration transition in a 
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proper way. Therefore, at present the discussion is restricted to so-called two-phase 
aproaches where the hadronic and the quark matter EoS are modeled separately 
followed by a subsequent phase transition construction to obtain a hybrid EoS. 

Widely employed for a description of quark matter in compact stars are ther- 
modynamical bag models of three-flavor quark matter with eventually even density- 
dependent bag pressure B(n), as in RefP^ 1 A qualitatively new feature of the phase 
structure appears in chiral quark models of the Nambu-Jona-Lasinio type which 
describe the dynamical chiral symmetry breaking of the QCD vacuum and its par- 
tial restoration in hot and dense matter, see RefpJ for a review. In contrast to 
bag models, in these approaches at low temperatures the light and strange quark 
degrees of freedom may appear sequentially with increasing densityJS-ES S o that 
strangeness may even not appear in the quark matter cores of hybrid stars, before 
their maximum mass is reached. Once chiral symmetry is restored, a rich spectrum 
of color superconducting quark matter phases may be realized in dense quark mat- 
ter, depending on it's flavor composition and isospin asymmetry® with far-reaching 
consequences for hybrid star phenomenology, e.g., M — R relationships and cooling 
behavior. 

We will consider here a color superconducting three-flavor NJL model with self- 
consistently determined density dependences of quark masses and scalar diquark 
gaps, developed in Refs.p^"® including the flavor-mixing KMT determinant inter- 
action JZDJ2Z} Only recently it became clear^'G^ that this flavor mixing is crucial 
for the possible stability of strange quark matter phases in hybrid stars. As a new 
aspect, we will couple the chiral quark sector to a Polyakov-loop potential which is 
a modification^ Q f the one used in standard PNJL modelsP^'^ 



§2. Modified PNJL model with color superconductivity 

We start from the path integral representation of the partition function for the 
modified color superconducting three-flavor PNJL model 

Z[T, V,{fl}] = j VqVq e ~ f ^m^d^-m-^ifl+iXa^q-C^+um ^-1) 

with the interaction Lagrangian £j nt = C qq + C qq , where 

8 

£qq = G S^2 [(<? r a<?) 2 + (q~H5T a q) 2 + G v (qij q) 2 

a=0 

- K [det/(g(l + 75)9) + det^l - l5 )q)] , (2-2) 
Cqq = G D y^ j (qij 5 T a X b Cq T )(q T Cij 5 T a X b q), (2-3) 

a,6=2,5,7 

with T a and Aj> being the antisymmetric Gell-Mann matrices acting in flavor and 
color space, respectively. We have suppressed the flavor index for the current quark 
mass matrix rh and the chemical potential matrix /}. The scalar diquark (Gd), 

vector (GV) and KMT (K) couplings are to be determined by hadron phenomenol- 
ogy, seeJ29 The Polyakov-loop <P = Tr c [exp(i/3A3(/)3)]/A^ c is an order parameter for 
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confinement, weighted with the phenomenological potential 

U(&) = (aT 4 + bfi 2 T 2 + C; u 4 )<2> 2 + a 2 T 4 ln(l - 6<2> 2 + 8<P 3 - 3<Z> 4 ) , (2-4) 

which is a modification 2 -' of the standard PNJL model, now accounting for an ex- 
plicit chemical potential dependence of U{<P), even at T = 0, which is not present in 
the traditional parametrizations by, e.g., RefsP^ 1 '^ 

Our description of quark matter is based on the grand canonical thermodynamic 
potential EDJIlHSDJliU which in the meanfield (stationary phase) approximation to 
the partition function (|2-ip is given by 



^ mf(T ' M) = 8GS + "T6G| SG~ V + AG~ D 



" / (0 E K + 2Tln { l + e ~ En/T )} - ^ + W (*o) , (2-5) 

whereE'n = E n {p, fi; (jlq, [i 3 , <f> u , <j> d , <j> s , u u ,uj d ,uj s , A ud , A us , A ds , <P ) are the quasi- 
particle dispersion relations, obtained by numerical diagonalization of the quark 
propagator matrix in color-, flavor-, Dirac- and Nambu-Gorkov spaces. The values 
of the meanfields (order parameters) are obtained from a minimization of f2{T, {/i}), 
which is equivalent to the selfconsistent solution of the set of corresponding gap 
equations. The subtraction of Qq garantees that in the vacuum i?(0, 0) = 0. For ap- 
plications to compact stars, we will have to add the contribution from leptons {e.g., 
electrons, muons and the corresponding neutrino flavors) and obey the constraints 
of color and electric neutrality as well as /3-equilibrium of weak interactions between 
quark flavors and leptons. 

The mean-field contribution of the Polyakov-loop potential could be viewed as a 
T, /^-dependent modification of the bag function AB{T,fx) = U{<Po) which accounts 
for possible changes in the pressure of the gluon sector related to a partial melting 
of the gluon condensate at finite T and fi. Therefore, it takes negative values. In 
the present contribution, we will discuss exploratory calculations of those effects 
for compact star structure, we will compare results with the above modified PNJL 
(mPNJL) model with those of approximating the influence of these effects by a 
simple (negative) bag constant AB. Note that the mPNJL model provides us with 
a more micoscopic picture on the origin of density (and temperature) dependence of 
nonperturbative QCD thermodynamics which in some phenomenological approaches 
are subsumed in a density dependent bag constant as, e.g., in RefA 2 ^ 

As pointed out in,^ due to the mixing of the light and strange flavor sectors by 
the KMT interaction, the difference in the critical chemical potentials for the chiral 
phase transitions in these sectors (which coincide with the onset of 2SC and CFL 
phases, respectively) gets diminished. This entails that the phase transition between 
hadronic matter (described by a realistic nuclear EoS, e.g., the DBHF one, see^ 1 ) 
and superconducting quark matter may eventually proceed directly into the CFL 
phase, provided the diquark coupling is sufficiently strong, see RefP^ 

In the present study, we will not investigate the dependence of the thermody- 
namics and phase structure on the strengths of the coupling constants which was 
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done in a previous workf^ 1 but rather adopt for their ratios with the scalar coupling 
strength the values given by the Fierz transformation of the one-gluon exchange in- 
teraction, i.e., Gy = 0.5 Gs and Gd = 0.75 Gs- The values of Gs, cutoff A and 
current quark mases m U) d and m s are chosen such as to describe the pion decay con- 
stant, pion and kaon masses as well as a light quark constituent mass of 350 MeV. 
The KMT coupling is fixed such that the rj — rj' mass splitting is obtained. 

In Fig. [1] we show the T = EoS as P((J>) for the ab initio DBHF approach to 
nuclear matter together with the results of evaluating the mPNJL model for quark 
matter in meanfield approximation (|2-5p for specific choices of the coefficient c in 
U(<I>o). It is interesting to note that those results compare well with those of the 
simplifying ansatz U{<£>o) = U{^q; \i = 0)+ AB, for the region close to the hadron-to- 
quark matter transition as relevant for applications to compact star structure. Note, 
that in the present model the range of Z\£?-values is limited to \AB\ < 50 MeV/fm 3 . 
For the mPNJL model, there is no such limitation of the values for the parameter c. 




Fig. 1. Equations of state for neutron star matter in beta-equilibrium. The hadronic phase is 
described by the DBHF EoS (solid line) and the quark matter EoS correspond to the mPNJL 
model (|2-5|l and its abridged version, for details see text. 



In Fig. [2] we show the M — R and M — n c relationships for nonrotating compact 
star sequences obtained as solutions of the Tolman-Oppenheimer-Volkoff equations 
for the hybrid EoS shown in Fig. [1] where the c-parameter of the mPNJL model (|2-4[) 
and the bag constant AB of the abridged mPNJL model are varied as respective 
free parameters. The diquark and the vector couplings are set to the values of the 
Fierz relation for one-gluon exchange, rjy = Gy/Gs = 0.5 and rjo = 0.75, resp. 
For Z\i3[MeV/fm 3 ] = —50 and —40 we obtain sequences of stable hybrid stars with 
2SC quark core, in the mass range 1.74 < M/M Q < 2.04 and 2.10 < M/M Q < 
2.12, respectively. For the other parameter choices, there is a direct transition from 
hadronic to CFL quark matter with a very small sequence of marginally stable 
CFL quark core hybrid stars at the respective maximum mass. A more detailed 
investigation is subject to ongoing work. 

The question arises whether flavor-mixing interactions other than the KMT one 
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have to be invoked. It seems necessary in a next step to include channels which 
appear after Fierz transformation of the KMT interaction and couple chiral conden- 
sates with diquark condensates.® This may lead to a new critical point in the QCD 
phase diagram® and, depending on the sign of the coupling, to a further reduction 
of the strange chiral condensate which enforces the flavor mixing effect studied here. 

The observational constraints for masses and radii are not yet settled. There is 
a lower limit for the M — R relation from RXJ 1856.5-3754® which requires either 
a large radius R > 14 km for a star with M = 1.4 M or a mass larger than ~ 2 M Q 
for a star with 12 km radius. The latter is accomodated with the present EoS. 
M — R relations for the quiescent binary neutron stars in globular clusters M13® 
and oj CeiF^* point to rather light and compact stars as described by those sequences 
obtained here, but without quark matter cores for masses M < 1.9 Mq. 

In order to test any of the predictions for quark matter phases in compact stars 
and their possible consequences, we look forward to future observational campaigns 
devoted to determine the masses® and M — R relationship for compact stars® with 
high precision and thus to constrain the dense matter EoSF^ 

§3. Conclusions 

The effect of density-dependent modifications of the gluon condensate on the 
sequential occurrence of superconducting two- and three-flavor quark matter phase 
transitions in the EoS for cold dense matter has been studied in a modified PNJL 
model with color superconductivity and flavor-mixing KMT determinant interac- 
tion. Different from previous work, here we fix the coupling constants from hadron 
phenomenology and Fierz relations which leaves the unknown chemical potential de- 
pendence of the modified Polyakov-loop potential as a free parameter of the present 
study. It is found that for sufficiently large values of the parameter c there can be a 
stable quark matter core in massive hybrid stars. 
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